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In this study, we have derived close form of the Casimir force for the non-interacting ideal Bose gas 
between two slabs in harmonic-optical lattice potential by using Ketterle and van Druten approx¬ 
imation. We find that Bose-Einstein condensation temperature T c is a critical point for different 
physical behavior of the Casimir force. We have shown that Casimir force of confined Bose gas 
in the presence of the harmonic-optical potential decays with inversely proportional to d s when 
T < T c . However, in the case of T > T c , it decays exponentially depends on separation d of the 
slabs. Additionally we have discussed temperature dependence of Casimir force and importance of 
the harmonic-optical lattice potential on quantum critical systems, quantum phase transition and 
nano-devices. 

PACS numbers: 05.30.Jp 


I. INTRODUCTION 

The Casimir effect is an attractive or repulsive long- 
range interaction between conducting boundaries [jj, 
which originates from the quantum mechanical fluctua¬ 
tions of the electromagnetic field in vacuum. The electro¬ 
magnetic field is confined in between plates with distance 
d, the fluctuations modes of electromagnetic fields need 
to have a node at the plate surfaces, so that only waves 
with wavelength A = nd/2 and integer n are permitted in 
the gap between plates. As a consequence the total field 
outside the gap produces a pressure on the plates which 
is higher than the one produced from inside, so the sur¬ 
faces are pushed together by force This effect has been 
experimentally measured in 1997 Q and 1998 [1] which 
confirms quantum field theory. It is shown that the force 
depends on size, temperature, geometry, surface rough¬ 
ness and electronic properties of the materials Q| (i . Its 
importance for practical applications is now becoming 
more widely appreciated in quantum field theory, gravita¬ 
tion and cosmology, Bose-Einstein condensation (BEC), 
atomic and molecular systems, mathematical physics and 
nano-technology. A good review on new developments in 
the Casimir effect can be found in Ref. Q. 

Now it is well-known that the fluctuations which can 
emerge from different physical mechanisms can lead to 
Casimir-like effect if it confined in the boundaries. For 
example, the thermal fluctuation in confined quantum 
critical systems can leads to Casimir-like effect. This 
emerges upon approaching a second-order phase transi¬ 
tion point which is characterized by the fact that due 
to the emerging collective behavior of the particles de¬ 
pends on thermal fluctuation. The effective force result¬ 
ing from the confinement of the fluctuations of the order 
parameter is called Casimir force. This phenomena firstly 
was predicted by Fisher and de Genes [7|. They showed 
that the fluctuation-induced force arise in critical systems 
which are close to the critical point, where dynamics are 
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governed by the thermal fluctuations of the order param¬ 
eter upon approach to the thermal phase transition. This 
effect later is called as the critical Casimir effect for quan¬ 
tum particle systems @]. Firstly in Ref. j&], the Casimir 
force between two slabs immersed in a perfect gas was 
calculated for various boundary conditions. It was found 
that the Casimir force has the standard asymptotic form 
with universal Casimir terms below the BEC critical tem¬ 
perature T c and vanishes exponentially above T c . After 
these seminal works, a great deal of effort has been de¬ 
voted to calculation of the Casimir force of the free and 
trapped Bose gas for different geometries and different 
boundary conditions flol - l35j . 

Another mechanism of Casimir-like effect arises from 
quantum size effect. It is expected that Casimir-like ef¬ 
fects may occur in geometrically confined quantum gas 
because of the wave character of the gas atoms at nano 
scale. Indeed it is shown that in the presence of bound¬ 
aries inside quantum particles can cause Casimir-like ef¬ 
fect if the bounded space is smaller than the maximum 
thermal de Broglie wavelength of the particles [3ll - l4^ |. 
Therefore, nowadays it has been attracted to Casimir 
force at nano scale because of the potential applications 
in nano-technology area since Casimir phenomena is a 
typically quantum size effect which appears at nano scale. 
In this area, researchers aim to produce many nano¬ 
device structures such as gas turbines, pumps, mixers, 
heat exchangers, valves, etc [38i - l4^ |. Therefore under¬ 
standing the Casimir effect on the thermodynamic be¬ 
havior of quantum gases confined within boundaries at 
nano scale in nano-devices is very significant problem. 
Because new nano-devices and technologies can be de¬ 
veloped based on the effects appearing at this scale. 

The subject of Casimir effect on thermodynamic be¬ 
haviors of gases at nano scale and in quantum critical 
systems is a new research area and may have many po¬ 
tential applications in nano technology and other area in 
physics. In this regard, discussing on the Casimir effect 
in the previous theoretical studies in Refs. (tUtI have po¬ 
tential to enlighten the thermodynamic behavior of con¬ 
fined quantum gases. So far Casimir force has been cal- 
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culated for either free or trapped Bose gas with harmonic 
potential loi-f30|. As it is known that the harmonic-optical 
potential which consists of combining harmonic and op¬ 
tical lattice potentials plays very important role in BEC 
systems as well as in other area of physics. The effects 
of harmonic-optical lattice potential on the BEC and its 
vortices and quantum phase transition have been consid¬ 
ered in Refs. [43l - [55j . But to our knowledge Casimir 
force of Bose gas trapped with harmonic-optical poten¬ 
tial has never been studied. In this study we investigate 
Casimir force of the Bose gas trapped with harmonic- 
optical potential between two infinite parallel slabs in 
the x — y direction that are separated by a distance d in 
the ^-direction. We obtain explicit form of the Casimir 
force of the ideal Bose gas trapped in harmonic-optical 
potential and we show that Casimir force decays with in¬ 
versely proportional to d 5 when T < T c while it decays 
exponentially for T > T c . 


II. MODEL AND ANALYTICAL RESULTS 


A. Theoretical framework 


In the original paper [lj Casimir force of vacuum 
fluctuation of electromagnetic field at zero temperature 
(T = 0) was defined as 


F c = --^[E(d)-E(oo) 


( 1 ) 


where E (d) is the ground-state energy (i.e. the vacuum 
energy) of the electromagnetic field in between the two 
conducting plates separated at a distance d. E (oo) is the 
energy at infinite. We consider the Casimir effect for the 
thermodynamical system of Bose gas between two infinite 
slabs. The geometry of the system on which some exter¬ 
nal boundary condition can be imposed is responsible for 
the Casimir effect. In recent studies it has been shown 
that the Casimir force for critical particle systems can 
be obtained from Casimir potential at finite temperature 
as 


B. Single particle energy 


In present model, we assume that ideal Bose gas is 
trapped in two dimensional harmonic-optical potential 
between two slabs separated at a distance d. If quantum 
particles move in three dimensional space in the absence 
external potential, single particle energy at finite temper¬ 
ature can be given by s (p x ,Py,Pz) = Px/^ m + Py/2m + 
p 2 /2m where m is the mass of the particle, and p x . p y and 
p z are the momentum along respectively x, y and 2 axis. 
However in our model particles are trapped in between 
quasi-two-dimensional slabs with two dimensional poten¬ 
tial. Hence the energy eigenvalues of particles at x and 
y-directions depend on the harmonic-optical potentials. 
For a simplicity, it can be assumed for present model 
that particle energy e depends on quantum number n XiV 
in x and y-directions and in the 2 -direction, energy of the 
particle is quantized as e nz = /(2md 2 ) where 

n z = 0,1,2,3,... and a indicates type of the boundary 
condition of the system. To define boundary conditions, 
it can be set as a = 1 and n z = 1,2,3,... for Dirich- 
let, a = 1 and n z = 0,1,2,... for Neumann, and a = 2 
and n z = 0, ±1, ±2,... for periodic. Therefore, the single 
particle energy for this model is arranged in the form 


— £rj 


ah 2 7 


2 md 2 


( 4 ) 


where e nx corresponds to two dimensional quantized 
energy of the particle in the x — y directions, cor¬ 

responds to energy of the particles for all degrees of free¬ 
dom. In order to obtain grand canonical potential of an 
trapped ideal Bose gas we have to compute single particle 
energy in Eq. (4) for the Bose gas. 

The first term of rhs in Eq. (4) must be determined to 
obtain exact expression for single particle energy e„ x a . 
The particle energy e n in x—y plane can be determined 
depending on harmonic-optical potential 


v (x, y) = V har (x, y) + Viat (x, y) (5) 


F c (T, p,d) =--^ [pc (T, p, d) - p c (T, p, 00 )] (2) 

where pc (T , p, d) and pc (T, p, 00 ) are the Casimir po¬ 
tential corresponds to fluctuations of thermally excited 
states of Bose gas at distance d and at infinite, respec¬ 
tively The Casimir potentials in Eq. (2) can be 

obtained from grand canonical potential of the particle 
system 

OO 

p (T, p, d) = -r 1 E In {1 ^ } (3) 

n 

where e n is the energy of Bose gas, which corresponds 
to energy eigenvalue of the particles. Additionally j3 = 
( kT ) 1 , fc is Boltzmann’s constant, T is temperature, and 
p is the chemical potential of the system. For simplicity 
we consider single particle energy instead of N particle 
Bose system. 


where harmonic potential and optical lattice potential 
are respectively given 

V h ar (x, y) = (uj 2 x 2 + u 2 y y 2 ) (6) 

Viat (x, y) = V 0 (sin 2 kx + sin 2 ky) (7) 

With m is the mass of the particle, are the fre¬ 

quencies of the harmonic potential along the coordinate 
directions x and y. Vo is the lattice potential depth and 
k = 2 tt/X' is the wave vector of the laser light, A' is the 
laser wavelength. In terms of k we can write the recoil 
energy e# = h 2 k 2 /2m {= huj^} as an energy scale for 
specifying the lattice depth. It is defined as the recoil 
energy that one atom requires when it absorbs one lat¬ 
tice photon [55| . 
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For the potential (5) it is impossible to find an exact 
analytical expression for the energy levels £ nx . An ap¬ 
proximate expression can be readily obtained by extend¬ 
ing the standard harmonic approximation for the opti¬ 
cal lattice potential [52|, [55]. Taylor expansion of Eq. (7) 
about the lattice site minimum at x = y = 0 leads to 

Vn 

Viat (x,y) = y{(l - COS 2 kx) + (1 - cos2fcj/)} (8) 
which is approximately equal to 

Viat (x,y) = V 0 k 2 {x 2 +y 2 ) - yfc 4 (a; 4 + y 4 ) . (9) 


As it can be seen from Eq. (9) that the first term can be 
organized in the form of the harmonic oscillator poten¬ 
tial \muj 2 at (x 2 +y 2 ) which yields an effective harmonic 
oscillator frequency 


Mlat 



zVErVq 

h 


( 10 ) 


with equivalent harmonic oscillator length l = 
sjh/mujiat■ la this case the localized states in the op¬ 
tical lattice can be approximated by a harmonic oscilla¬ 
tor states, hxiat . Under this consideration the effective 
confining potential can be given by 


Veff {x, V) = \rn [(ul + wf at ) x 2 + 


This potential is characterized by a single particle energy 
levels with Dirichlet boundary condition given by 

^nx,y - ’H'xfyXx A T A^yj (12) 

where w' = y/u>f + cof at with i stands for x and y, on the 
other hand, n x , n y = 0,1, 2,... and UJi = | (u> x + ui y ) is 
the mean of combined frequencies. In deriving Eq. (12), 
the second term in Eq. (11), quartic term is treated per- 
turbativly by using the normal ladder operators (52| . 
This term provides a shift in the oscillator state ener¬ 
gies. Up to the first order term this shift is given by 

[(n 2 + n 2 y ) + (n x + n y ) + l] (13) 

where u>r = hk 2 /2m is the recoil frequency. Note that 
for a fixed value of er = Hlor the quantity Ae n grows 
rapidly with n' s. This is not surprising, since the higher 
excited states of the harmonic oscillator extend to larger 
and larger values of x and y. In such case the pertur¬ 
bation term can have important effect for higher excited 


K + “lat) v 2 ] - Y‘ fc4 (* 4 + y4 "> ■ ( n ) 


state. Also from Eq. (13) we see that the higher the value 
of n x and n y are smaller than the values of er for which 
reliable values may be obtained from the experimental 
set-up conditions, er = h 2 k 2 /2m, where k is determined 
from the laser beam wavelength. Thus in order to keep 
the value of £r in its experimental range one has to ne¬ 
glect the term which contains (n 2 + n 2 ) in Eq. (13). Un¬ 
der this approximation the single particle energy level 
depends on n x and n y is given 

= n x H (u' x - y) + n y H (w y - y) + e ° ( 14 ) 

where Eq = ktix is the single particle ground state energy 
with ol = TjJ\ — is the meaning of the combined frequen¬ 
cies. Eq. (14) means that particle oscillates in quasi-two 
dimensional space. On the other hand, there is an energy 
contribution to particle in the 2 -direction. For Dirich¬ 
let boundary condition this contribution can be given as 
M; n 2 . Hence the single particle energy in all directions 
can be defined by 


= n x h I 


r* ~ T) +nyh K “ + w + 


£o 


(15) 


Now with help Eq. (15) grand canonical potential for two dimensional trapped Bose gas in harmonic-optical potential 
in between two slabs can be obtained. 


C. Thermodynamical Potential 

Grand canonical potential (3) can be written as 

OO 

p ( T , )jl, d) = —ipo — kT In (l — ze~P En *,y,z ) (16) 

n=1 

where first term ip q is ground state potential and second term corresponds to the potential of excited particles. We 
have neglected the ground-state term which does not give any contribution to Casimir potential at finite temperature 
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since the Casimir effect was caused by thermal or near-critical fluctuations of the order parameter upon approaching 
the condensation transition. The second term in Eq. (16) for finite particle can be presented as 

OO OO j 

lp (T,v,d) = -kTj2Y,- e ~ jllSnx ’ v ” U?) 

n= 1 j=l ^ 

If Eq. (15) is put into Eq. (17), grand canonical potential for two dimensional trapped Bose gas can be written in 
generalized form 


¥> (T , n, d ) 


OO OO OO OO j s / 

- kT EE E E 7 ex p {\ n * h (E 

n x n y n 2 = lj=l J X 




+ 


2 md 2 


(18) 


where z = e Here we have neglected the energy £q in Eq. (15) since ground state does not contribute to Casimir 
potential. Eq. (18) can be simplified in the form 


<P (T, n, d) 


-XT E E 


' exp { —jP (EcF n z) } 


=i ~r[ I 1 - ex P {iP h K - }) (! - ex P {jP* K - }) 


(19) 


To evaluate the sum over j, we need an approximation due to Ketterle and van Druten [56[ which yields near 
exact results. This approximation uses the thermodynamic limit fax <C kT to write (l — exp {j/3H (w( — ^)}) = 
j(3H (oj( — ^) where i = x,y. The sum over j has the factor 1 /j 2 and is convergent. Thus, Eq. (19) is written as 


<p( T ,f* 1 d) = ~kT- 


kT \ “ 
h ) 


K-¥) 


EE^exp \-jP 


\2md? 


( 20 ) 


where we set n z = n for simplicity. This analytical expression corresponds to the grand canonical potential of trapped 
Bose gas in quasi-two-dimensional harmonic-optical potential in between two slabs due to Ketterle and van Druten 
approximation. 


D. Casimir potential 

Casimir potential can be obtained from Eq. (20). How¬ 
ever, Eq. (20) includes contributions of the bulk, surface 
and fluctuation of excited states of the Bose gas between 
two slabs separated by d. Bulk and surface potentials 
do not contribute to the Casimir force, because the force 
due to the bulk term is counterbalanced by the same 
contribution acting from outside the slabs when they are 
immersed in the critical medium 1,0,113 , on the other 
hand, the surface potential caused from two slab-boson 
gas interfaces, does not change with changing d. However 


thermal fluctuations of excited states of Bose gas leads 
to Casimir potential. Therefore, to obtain Casimir po¬ 
tential, these contributions must be se par ated into com¬ 
ponents. By using the Jacobi identity [53 



where b = jn(X/d) 2 /2 with b > 0 and A = 
is thermal de Broglie wavelength of the particles, these 
contributions can be separated as 


f {T, d) = - 


(kT) 3 y Z l 

v K - W) 


2 Vb 


+ 


Vb 


E< 

71=1 


7r n 2 /b 


( 22 ) 


In Eq. (22) the first term corresponds to bulk potential fbuik (T,f,d) and second term is the surface potential 
fsurf (T,f,d). However, the last term in Eq. (22) corresponds to the Casimir potential for combined potential 

\3 


fc (T , n,d) = - 


(kry 


v K - ¥) K-^VttA 


2 d “ “ e j0T 


EE 

1=1 n =1 


7/2 


-2(nd/\) /j 


(23) 


where sum over j and n. The summation in right hand side of Eq. (23) can be converted to integral in the limit 

d/i < 1, 


EE 


eiPk 

JJ2 


e -2(nd/X) 2 /j 



v -6 e -px 2 - q lx 2 dx 


(24) 
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where p = u 2 /2 = (—2/3/j.) 1//2 d/X ~ d/£, £ is the correlation length, q = 2n 2 and a; 2 = (-j|) 2 r. Thus, Casimir potential 


is 


{T,p,d) = -- 


2 (fcT) 3 


2 /A 


V (fa~^) K - if) V 7T Vd 

After some mathematical algebra, Casimir potential can be obtained as follows 

3 ( kT ) 3 


-6 e - P x 2 -qlx 2 dx 


ipc (T, p,d) = - 


64 h 2 i: 2 fa x -^) fa-**) \d 


\\ 4 /I _L 0».», _L /I.,2^2 


E 


1 + 2un + 4 u 2 ?z 2 /3^ e _ 2u „ 


If A = hy/P/m is inserted into Eq. (26), Casimir potential is given in a simple form 


pc {T,p, d) = -- 


3H 2 kT 


E 


1 + 2wn + 4u 2 n 2 /3^ ^ 2u „ 


64tt 2 to 2 fa - **) fa y - **) d 4 
which originates from fluctuations of thermally excited states of Bose gas in between slabs. 


(25) 


(26) 


(27) 


E. Casimir Force 


J 


Now we can compute Casimir force for tapped ideal Bose gas in harmonic-optical potential between two slabs by 
using of Eqs. (2) and (27). Casimir potential in Eq. (27) goes to zero i.e., pc (T, p, oo) —> 0 for d —> oo. Therefore 
Casimir force in Eq. (2) reduces to simple form 

F c (T,p,d) = --^p c (T,p,d) (28) 

which implies that Casimir force of Bose gas changes due to the changing of distance between slabs. From Casimir 
potential (27), an general analytical form for the force is given simply by 


F c (T, p, d) 


3kTH 2 


\ ' 


16tt 2 to 2 fa x - fay - ^t) d 5 


1 + 2 un + 4:U 2 n 2 /3\ un 


(29) 


However this expression does draw a physical picture of 
the Casimir force for the Bose gas since it depends on 
chemical potential p of the gas. It is well known that p 
has a curious behavior depends on temperature. In the 
Bose gas, thermal fluctuations cause second order phase 
transition which is called BEC occur at critical tempera¬ 
ture T c . Chemical potential p equal to zero at the T = T c 
and p < 0 when T > T c . It is known that correlation 
length £ of the thermal fluctuations diverge to infinite at 
the second order phase transition temperature T — T c 
where most of the particles begin to collapse to ground 
state and occur condensation. However contribution to 
Casimir potential does not come from ground state, it is 
meaning that BEC does not cause Casimir effect. The 
origin of the Casimir effect in the critical particle systems 
is excited states. The correlation length £ of thermal fluc¬ 
tuation of the excited states of the quantum gas is finite 
when T gl T c . If the geometry is restricted to be that 
of two parallel plates separated by distance d in the , 2 - 
direction to be k = irn/d (for the Dirichlet boundary 
condition), correlated fluctuations become dominate and 


the free energy of the Bose gas is altered in the medi¬ 
ate vicinity of the critical point d/£ < 1 leading to the 
Casimir force. Therefore, it is suggested that Casimir ef¬ 
fect for critical particle systems is caused by thermal or 
near-critical fluctuations of the order parameter upon ap¬ 
proaching the condensation transition. In this picture, we 
assume that in the near T gj T c , the chemical potential to 
be p = 0. In this case, when T ~ T c , z = 1 behaves as the 
bulk critical point for Casimir force. Hence for p = 0, the 
sum in Eq. (29) takes £ (5) = Therefore, the 

Casimir force in the vicinity the BEC transition tempera¬ 
ture T c is obtained in terms of Rieamann’s Zeta function 
as follow 


F c (T, p, d) = 


3kTh 2 ((5) 30) 

167T 2 m 2 fa x - **) fa y - **) (r 


where fa = ^fa+fa at and fa = yjfa+fa at - For 
isotropic case, the Casimir force in combined potential 
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reduce to 

F c (T, n, d) 


3kTh 1 2 C (5) 1 

167t 2 4 to 2 ( u )' — ^) 2 


(31) 


where w' = ^/co 2 + uj‘f at . It can be seen that Casimir force 
in BEC phase or T ~ T c decays inversely proportional to 
d 5 6 7 8 9 10 for harmonic-optical potential as well as in the case 
of harmonic potential. However, as it can be seen from 
Eqs. (30) and also (31) Casimir force of Bose gas clearly 
depends on frequencies of harmonic and optical lattice 
potentials. 

On the other hand, in the case of T > T c , the chemical 
potential is n ^ 0 and the double sums in Eq. (23) for 
d/£ 3> 1 satisfy the condition 


E 


j =i 


e /3/b 

w 


E 


e -2(ud/\) 2 /j 


n= 1 


< 


C(7/2) _ 

gs/—8fid/\ _ ^ 


(32) 


Therefore in the limit of large d , Casimir potential for 
harmonic-optical potential can be given 

<pc(T,Li,d)~e-^ d / x . (33) 

As it can be expected that for d/£ 1 the Casimir po¬ 

tential has exponential form which indicates that Casimir 
force also decays exponentially with d. Additionally, at 
the high-temperature limit; the correlated fluctuations 
between the two slabs are no longer long-ranged there¬ 
fore the confi ning boundaries are subject to a vanishing 
Casimir force [Tot fl7|. 


III. CONCLUSION 

In this study, we have derived close form of the 
Casimir force for the ideal Bose gas between two slabs 


in harmonic-optical potential using by Ketterle and van 
Druten approximation. We have shown that Casimir 
force of confined Bose gas in the presence of the 
harmonic-optical potential decays with inversely propor¬ 
tional to d 5 when T <T C . However, in the case of T > T c , 
it decays exponentially depending on separation d of the 
slabs. Analytical results indicate that BEC critical tem¬ 
perature T c plays important role on the defining of the 
Casimir force in confined critical particle systems. As dis¬ 
cussed above, when fluctuations are long ranges (compa¬ 
rable to smallest dimension of the system) near the BEC 
transition in Bose gas there will be a Casimir force be¬ 
tween the confining surfaces. Based on above discussing 
we can suggest that similar behavior may happen for any 
confined system near a critical point if the order param¬ 
eter has critical fluctuations on a scale £ which becomes 
large, diverging in the bulk/thermodynamic limit. Fi¬ 
nally, we state that the role of the Casimir effect in con¬ 
fined systems such as BEC, quantum phase transition at 
ultra cold Bose systems and nano structures is still open 
problems. Harmonic-optical lattice potential has a sig¬ 
nificant role on these systems. We hope that the present 
analytical results may contribute to theoretical and ex¬ 
perimental studies in these area. 
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